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ABSTRACT 

We derive the growth of SMBHs relative to the stellar content of their host galaxy predicted 
under the assumption of BH accretion triggered by galaxy encounters occurring during their 
merging histories. The latter are described through Monte Carlo reaUzations, and are con- 
nected to gas processes, star formation and BH accretion using a semi-analytic model of 
galaxy formation in a cosmological framework. This allows us to connect the star formation 
process in the host galaxies to the growth of Supermassive Black Holes. We show that, within 
this framework, the ratio T = {MBH/Mt)iz)/iMBH/M,){z = 0) between the Black Hole mass 
and the galactic stellar mass (normaUzed to the local value) depends on both BH mass and 
redshift. While the average value and the spread of r(z) increase with z, such an effect is larger 
for massive BHs, reaching values F 5 for massive Black Holes (M > lO'^ Mq) at z > 4, 
in agreement with recent observations of high-redshift QSOs; this is due to the the effective- 
ness of interactions in triggering BH accretion in high-density environments (where massive 
haloes form) at high redshifts. 

To investigate how different observations of r(z) fit within our framework, we worked out 
specific predictions for sub-samples of the simulated galaxies corresponding to the different 
observational samples for which measurements of F have been obtained. We found that for 
Broad Line AGNs at intermediate redshifts 1 < z ^ 2 values of F « 2 are expected, with a mild 
trend toward larger value for increasing BH mass. Instead, when we select from our Monte 
Carlo simulations only extremely gas rich, rapidly star forming galaxies at the epoch of peak 
in the cosmic star formation (2 < z < 3), we find low values 0.3 < F < 1.5, consistent with 
recent observational findings on samples of sub-mm galaxies; in the framework of our model, 
these objects end up at z = in low-to-intermediate mass BHs (M < 10^ Mq), and they do not 
represent typical paths leading to local massive galaxies. The latter have formed preferentially 
through paths (in the M. - Mbh plane) passing above the local M* - Mbh relation. We discuss 
how the global picture emerging from the model is consistent with a downsizing scenario, 
where massive BHs accrete a larger fraction of their final mass at high redshifts z > 4 
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1 INTRODUCTION 

The discovery of tight local correlations between the mass Mbh of 
Supermassive Black Holes (SMBHs) and global properties of their 
host galaxies, like the stellar mass M, or the velocity dispersion cr,, 
(Kormendy & Richstone 1995; Magorrian et al. 1998; Ho 1999; 
Gebhardt et al. 2000; Ferrarese & Merritt 2000; Marconi & Hunt 
2003; Haring & Rix 2004; Kormendy & Bender 2009) constitutes 
an important breakthrough in the understanding of galaxy evolu- 
tion. If the growth of SMBHs is considerably contributed by accre- 
tion over cosmological times (Soltan 1982; Yu & Tremaine 2002; 
Marconi et al. 2004; Merloni & Heinz 2008) this implies that most, 
if not all, galaxies have hosted an Active Galactic Nucleus (AGN) 
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in the past, and that a strong physical connection exists between 
galaxy formation and the growth of SMBHs. 

Understanding the mechanisms responsible for such a con- 
nection, and establishing the relative time scales for star forma- 
tion/assembly and for SMBHs growth, requires the measurement 
of the above correlations at higher redshifts. For such a scope, the 
most straightforward correlation to study is the Mgjj - M, relation, 
since the other fundamental relation Mgg - cr, would require the 
measurement of stellar velocity dispersions which are very difficult 
to achieve at z > 1 in normal galaxies. 

The Mbh - M, relation in the local Universe has been widely 
investigated in the literature. Magorrian et al. (1998) were the first 
to find a correlation between the black hole mass and the stellar 
mass, however they found a very large scatter (~0.51 dex). This 
relation was then re-examined using more reUable Mbh and M, 
measurements (Merritt & Ferrarese 2001, Marconi & Hunt 2003, 
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Haring & Rix 2004); these studies showed that the scatter in the 
^BH ~ relation is comparable to the scatter in the Mgu - cr, 
relation {~0.3 dex, see Haring & Rix 2004). It is also worth noting 
that the relation is BH mass dependent (Haring & Rix 2004), 
although the mean value of the Mbh/M, ratio obtained from the 
Haring & Rix (2004) sample: {tog(MBH/M,)} ^-2.8, is consistent 
with those derived by Merritt & Ferrarese (2001) and Marconi & 
Hunt (2003) when the latter is lowered by a factor of 5/3 to account 
for more recent estimates of the bulge virial masses (see Maiolino 
et al. 2007). 

The SMBH masses at z > 0.5 are generally estimated in lu- 
minous type 1 AGNs like quasars (QSOs) and radio galaxies; for 
example, virial BH mass estimators are based on the assumption 
that the broad-line region of the AGNs are dominated by the grav- 
ity of the SMBH (see, e.g., Peterson & Wandel 2000), and provide 
mass estimates of SMBHs at higher redshifts (see Wandel, Peterson 
& Malkan 1999; Kaspi et al. 2000; McLure & Dunlop 2002, Mer- 
loni et al. 2009), extending out to z ~ 6 (e.g., Willott, McLure & 
Jarvis 2003; McLure & Dunlop 2004; Vestergaard et al. 2004; Wal- 
ter et al. 2004; Riechers et al. 2008). The evolution of the Mbh/M, 
ratio from the local value is generally quantified in terms of the 
parameter; 

r = {Mbh/M,){z)/{Mb„/M,){z = 0). (1) 

Although the AGN activity makes it difficult to derive stellar 
masses for the host galaxies at z > 1.5, several observational 
works suggest that at such redshifts the mean Mbh/M, ratio is 
significantly higher than the present value, corresponding to F > 1 . 

Peng et al. (2006), applying the virial method on a sample 
of 31 gravitationally lensed AGNs and 20 non-lensed AGN and 
deriving the host galaxy stellar mass from the R-band luminosity, 
found the Mbh/M, ratio to be higher than the present value by a 
factor F « 2 for AGNs at 1 < z < 1.7 and by a factor of F ~ 4 at z > 
1.7; analysing the hosts of 89 broad-line AGNs in the COSMOS 
survey Merloni et al. (2009) derive a positive evolution fitted with 
F K (1 + z)"^'' in the redshift range 1 < z < 2.2; exploiting radio- 
loud unification McLure et al. (2006) obtained virial (linewidth) 
BH mass estimates from the 3C RR quasars, and the stellar mass 
estimates from the 3C RR radio galaxies, thereby providing black 
hole and stellar mass estimates for a single population of early-type 
galaxies, finding F « (1 + z)^ for z ;S 2; at even higher redshifts 
z > 4, measurements of the dynamical mass from the CO emission 
of molecular gas in the hosts of high-redshift QSOs (Walter et al. 
2004; Riechers et al. 2008) yield large values of F = 5 - 10. 

While the above studies indicate that the growth of SMBHs 
is faster than the stellar mass assembly, they all focus on lumi- 
nous AGNs, and are thus biased towards selecting the most mas- 
sive SMBHs (see Lauer et al. 2007). In contrast to the above re- 
sults, from the analysis of a sample of sub-mm selected galaxies 
(SMGs) Borys et al. (2005) find F to decrease with redshift for 
z < 2; more recently, Alexander et al. (2008) found F(z = 2) = 0.3 
for SMGs at z ~ 2 , indicating that for such objects the growth of 
SMBHs actually lags that of the host stellar mass. These latter stud- 
ies are affected by a different bias, in that they select ultra-luminous 
(L ^ 10'^ Lg), gas-rich galaxies (with gas fraction relative to bary- 
onic mass fgas ^ 0.6, sec Grcvc et al. 2005; Tacconi et al. 2008), 
characterized by the conversion of a significant fraction of their ini- 
tial gas reservoir of 10'" - 10'' Mq into stars in a short times scale 
of a few 10* yrs (see Tacconi 2008 and references therein); given 



the correlation between star formation rate and stellar mass, such 
observations tend to be biased towards massive stellar hosts. 

Attempts to derive the evolution of the global <F(z)) averaged 
over for whole AGN population requires a different approach, based 
on integrated observables. By comparing the redshift evolution of 
the integrated BH and of the stellar mass densities Merloni, Rud- 
nick & Di Matteo (2004) found that (F) « (1 + z)" with or ^0.5. 
Such a results is consistent with the upper limits on the evolution 
of (F) derived by Hopkins et al. (2006) for z g 2. 

Thus, a major challenge for next studies of AGN and galaxy 
formation is understanding which of the above classes of paths F(z) 
is the dominant one in the cosmological evolution of galaxies and 
SMBHs, and what are the physical mechanisms determining the 
particular path followed by a galaxy as a function of its properties 
in different cosmic epochs. While on the observational side such 
a task would require a large, unbiased sample of AGNs, a com- 
plementary approach is to consider cosmological models of galaxy 
and SMBH evolution which assume a given physical mechanism 
for the AGN feeding and for its connections with the evolution of 
the host galaxy; then, comparing the predictions of such models 
for specific classes of objects with different existing observations 
would constitute a powerful probe for the proposed physical link 
between galaxies and SMBHs. 

Indeed, in the recent years, significants steps forward had been 
taken toward the modelling of the cosmological co-evolution of 
galaxies and AGNs, starting from the seminal papers by Silk & 
Rees (1998); Fabian (1999); Cavaliere & Vittorini (2000); Wyithe 
& Loeb (2003). On the one hand, recent high-resolution hydrody- 
namic N-body simulations (see, e.g., Di Matteo, Springel & Hem- 
quist 2005, Springel 2005; Hopkins 2005, 2006) have shown the 
importance of galaxy mergers as triggers for AGN accretion, and 
the role of the AGN energy feedback in the subsequent evolution 
of the host galaxy. However, such simulations are affected by a 
limited exploration of the sub-grid prescriptions for the physics of 
SMBHs; most important, they necessarely focus on specific galaxy 
systems, making it difficult to assess the statistical relevance of a 
particular path in the building up of the local Mbh - M, relation. 
On the other hand, semi-analytic models (SAMs, Monaco, Salucci 
& Danese 2000, Kauffmann & Haenhelt 2000; Volonteri, Haardt & 
Madau 2003; Granato et al. 2004; Menci et al. 2006; Croton et al. 
2006; Bower et al. 2006; MaruUi et al. 2008) can probe the whole 
set of possible paffis followed by galaxies in the Mbh - M, plane 
during their evolution. Since their results depend on the mecha- 
nisms assumed for the trigger of the AGN accretion and for the 
AGN energy feedback, SAMs are suitable to probe the statistical 
outcomes provided by the implemented physical processes Unking 
SMBHs with ffieir host galaxy. 

Here we adopt a state-of-the-art SAM (Menci et al. 2006, 
2008) which implements a specific, physical model for the BH 
accretion based on interactions as triggers for both the fueling of 
SMBHs and the star bursts. In turn, the interaction rate, the amount 
of gas available for BH accretion, and the stellar content of galax- 
ies are computed from Monte Carlo realizations of galaxy merg- 
ing trees, enabling us to follow the evolutionary paths of galaxies 
and of their SMBHs over cosmological times, from the collapse of 
their progenitors from the primordial density field to the present. 
Based on such a model we derive the evolutionary paths for the 
growth of SMBHs relative to that of ffie stellar mass of their host 
galaxies, to investigate which are the dominant paths followed by 
different galaxy populations in building up the local Mbh - M, re- 
lation in a scenario characterized by AGN fueling driven by galaxy 
interactions. We then compare different model predictions for the 
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MgnlMt ratio of luminous intermediate- and high-redshift galax- 
ies, and for intermediate redshift SMGs witli existing observations, 
and discuss the physical origin of the different paths T(z) followed 
by galaxies and SMBHs during their evolution. Finally, we shall 
provide specific predictions typical of the scenario characterized by 
interaction-driven AGN fueling, to be tested with next experiments. 



2 THE MODEL 

We use the semi-analytic model as is described in details in Menci 
et al. (2005, 2006, 2008); this connects, within a cosmological 
framework, the accretion onto SMBHs and the ensuing AGN activ- 
ities with the evolution of galaxies. Here we recall the basic points. 



2.1 Hierarchical Galaxy Evolution 

Galaxy formation and evolution is driven by the collapse and 
growth of dark matter (DM) haloes, which originate by gravita- 
tional instability of overdense regions in the primordial DM den- 
sity field. This is taken to be a random, Gaussian density field with 
Cold Dark Matter (CDM) power spectrum within the "concordance 
cosmology" (Spergel et al. 2007) for which we adopt round param- 
eters fl/i = 0.7, Q.^) = 0.3, baryonic density Q./, = 0.04 and Hubble 
constant (in units of 100 km/s/Mpc) h = 0.7. The normalization of 
the spectnun is taken to be erg = 0.9 in terms of the variance of the 
field smoothed over regions of 8 Mpc. 

The merging rates of the DM haloes are provided by the Ex- 
tended Press & Schechter formalism (see Bond et al. 1991; Lacey 
& Cole 1993). The clumps included into larger DM haloes may 
survive as satellites, or merge to form larger galaxies due to binary 
aggregations, or coalesce into the central dominant galaxy due to 
dynamical friction; these processes take place over time scales that 
grow longer over cosmic time, so the number of satellite galaxies 
increases as the DM host haloes grow from groups to clusters (see 
Menci et al. 2005, 2006). 

The processes connecting the baryonic components to the 
growing DM haloes (with mass m) are described in our previous 
papers (e.g., Menci et al. 2005). The gas at virial equilibrium with 
the DM potential wells undergoes radiative cooling. The cooled gas 
mass rric settles into a rotationally supported disc with radius 
(typically ranging from 1 to 5 kpc), rotation velocity and dynam- 
ical time tj = r^/vj. Two channels of star formation may convert 
part of such a gas into stars: 

i) quiescent star formation, corresponding to the gradual conversion 
into stars at a rate m. oc nic/td 

ii) starbursts triggered by interactions, which destabilize part of the 
cold gas available by inducing loss of angular momentum. Note 
that galaxy interactions (among galaxies with relative velocity Vrei) 
include not only merging but also fly-by events. Part of energy re- 
leased by SNae following star formation is fed back onto the galac- 
tic gas, thus returning part of the cooled gas to the hot gas phase 



2.2 Accretion onto SMBHs and AGN emission 

The model also includes a treatment of SMBHs growing at the cen- 
tre of galaxies by interaction-triggered inflow of cold gas, foUow- 
ing the physical model proposed by Cavaliere & Vittorini (2000) 
and implemented in Menci et al. (2006, 2008). The accretion of 
cold gas is triggered by galaxy encounters (both of fly-by and of 
merging kind), which destabiUze part of the cold gas available by 



inducing loss af angular momentum. The rate of such interactions 
is 



: riT £(r,, V, Vrel) Vrel- 



(2) 



Here iij is the number density of galaxies in the same halo, is 
their relative velocity, and 2 is the cross section for such encoun- 
ters which is given by Saslaw (1985) in terms of the tidal radius r, 
associated to a galaxy with given circular velocity v (see Menci et 
al. 2003, 2004). 

The fraction of cold gas accreted by the BH in an interaction 
event is computed in terms the variation A j of the specific angular 
momentum j x Gm/vj of the gas to read (Menci et al. 2003) 
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Here h is the impact parameter, evaluated as the average distance of 
the galaxies in the halo. Also, m' is the mass of the partner galaxy 
in the interaction, and the average runs over the probability of find- 
ing such a galaxy in the same halo where the galaxy with mass m is 
located. The values of the quantities involved in the average yield 
values of face :£ 10"^. For minor merging events and for the en- 
counters among galaxies with very unequal mass ratios m' <g: m, 
dominating the statistics in all hierarchical models of galaxy for- 
mation, the accreted fraction takes values 10"' < face % 10"^. 

The average amount of cold gas accreted during an accretion 
episode is thus Ainacc = face '"o and the duration of an accretion 
episode, i.e., the timescale for the QSO or AGN to shine, is as- 
sumed to be the crossing time t = rj/v^ for the destabilized cold 
gas component. 

The time-averaged bolometric luminosity so produced by a 
QSO hosted in a given galaxy is then provided by 



(4) 



We adopt an energy-conversion efficiency rj - 0.1 (see Yu & 
Tremaine 2002), and derive the X-ray luminosities Lx in the 2- 
10 keV band from the bolometric corrections given in Marconi 
et al. (2004). The SMBH mass msH grows mainly through accre- 
tion episodes as described above, besides coalescence with other 
SMBHs during galaxy merging. As initial condition, we assume 
small seed BHs of mass 10^ Mq (Madau & Rees 2001) to be ini- 
tially present in all galaxy progenitors; our results are insensitive to 
the specific value as long as it is smaller than some 10' Mq. 



2.3 AGN feedback 

Finally, our SAM model includes a detailed treatment of AGN feed- 
back. This is assumed to stem from the fast winds with velocity up 
to lO^'c observed in the central regions of AGNs (Weymann, Car- 
swell & Smith 1981; Turnshek et al. 1988; Crenshaw et al. 2003; 
Chartas et al. 2002; Pounds et al. 2003, 2006; Risaliti et al. 2005); 
these are usually though to originate from the acceleration of disc 
outflows due to the AGN radiation field (Proga 2007 and references 
therein, Begelman 2004). These supersonic outflows compress the 
gas into a blast wave terminated by a leading shock front, which 
moves outwards with a lower but still supersonic speed and sweeps 
out the surrounding medium. Eventually, this medium is expelled 
from the galaxy. 

Quantitatively, the energy injected into the galactic gas in such 
inner regions is taken to be proportional to the energy radiated by 
the AGN, AE = €agn 1 <^ Antacc- The value of the energy feedback 
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Figure 1. Left Panel: The local predicted Mbh - relation is compared with data by Haring and Rix (2004, diamonds), and Marconi & Hunt (2003, squares, 
here M, is derived using the best-fitting virial relation of Cappellari et al. 2006); the colour code represents the density (per Mpc^) of BHs in a given Mbh ~ M, 
bin, as indicated by the upper colour bar. The contour plot allows to show how the predicted behaviour of the scatter (growing with decreasing stellar mass) is 
consistent with the observed scatter of the data points. 

Central Panel: The predicted Mbh ~ M, relation at z = 4 (colour coded as above) is compared with data corresponding to observations of broad-line AGNs 
and host galaxies in the range 4 < z < 6; details on the data and corresponding references are provided in the Appendix A. We also show as a solid line the 
predicted median local value, corresponding to the median local value of Mbh for model galaxies with given M» . 

Right Panel: We show some of the paths in the Mg//(/) - M,(t) plane followed, during their evolution, by BHs (and by their host galaxies) reaching a final 
mass of Mbh(z = 0) > 10'° Mq. 



efficiency for coupling with the surrounding gas is taken as exojv = 
5 IQ-^ (see Menci et al. 2008). 



3 RESULTS: OVERVIEW OF SMBH GROWTH 

As a first step in the study of the relative growth of BHs and host 
galaxies, we show in fig. 1 the local (left panel) and the high- 
redshift (z = 4, central panel) Mbh - M, relation that we obtain 
from our model. 

The predicted local relation is consistent with observations, 
and the data lie on the predicted confidence region represented by 
the contour plot, although the predicted distribution shows a small 
offset from the observed best-fitting relations derived by Marconi 
& Hunt (2003) and Haring & Rix (2004). Although, in principle, 
the model parameters (like the AGN feedback efficiency e^Qf^ or 
the normalization of the star formation efficiency) could be tuned 
to optimize the fitting to the observed local Mbh - ^* relation, 
this would also affect the predicted properties of the galaxy pop- 
ulation (e.g., the evolution of the luminosity function, the colour 
distributions, the local TuUy-Fisher relation) at both low and high 
redshifts. Since the model is intended to provide a unified descrip- 
tion of AGNs and galaxy evolution, we chose to adopt the same 
fiducial model adopted in our previous works since it provides a 
good match to a wide set of galactic properties (see Menci et al. 
2006, 2008). Instead, a robust feature of the Mbh - M, is consti- 
tuted by the behaviour of the scatter, which is almost independent of 
the model parameters and even on the specific mechanism assumed 
to trigger the AGN feeding. In fact, its decrease with increasing 
Mbh 's typical of hierarchical scenarios, and it is due to the early 
assembly of progenitors into a unique main progenitor which char- 



acterizes the merging histories of massive objects formed in biased, 
high density environments (see Menci et al. 2008). 

Note that the Mbh - M» relation is predicted to evolve with z 
as shown by the central panel of fig. 1, with the massive end of the 
distribution reaching the region corresponding to the observations 
of luminous, broad-line AGNs at z > 4 (see Appendix A for details 
on the data). A typical prediction of our interaction-driven model is 
that the evolution of the Mbh - M, relation increases for increasing 
stellar or BH mass; indeed, while low-mass BHs in model galaxies 
at high redshifts are characterized by Mbh/M, ratios close to their 
local values (at least for the majority of them), for the most massive 
simulated BHs the Mbh/M, is substantially larger than the model 
local values. We can quantify the evolution of the Mbh/M, in our 
model through the quantity T defined in Sect. 1 (eq. 1); since our 
model predicts a whole distribution of local Mbh for each stellar 
mass M, (not just a simple power-law, as is generally adopted to 
fit the local data), we compute T as the deviation of the Mbh/M, 
ratio of each model black hole from its own local value. Then for 
galactic hosts with stellar mass M, « 10" Mq the BH masses at 
z > 4 are typically Mbh ~ 10'"' Mq (see central panel of fig. 1), 
more massive than the predicted local value by a factor F a; 5. These 
objects constitute the progenitors of the local extremely massive 
BHs, and their growth with time is represented by the paths shown 
in the right panel of fig. 1 . Note how such paths are characterized 
by an assembly of BH masses extremely rapid in the early phases, 
which is faster than the stellar mass growth since they approach 
their final position in the Mbh - ^» plane passing above the local 
Mbh - M, relation. 

The physical origin of the above behaviour can be understood 
as follows; such high-redshift massive BHs (the counterpart of 
those observed by Walter et al. 2004; Maiolino et al. 2007; Riech- 
ers et al. 2008) are formed from galaxies collapsed in biased, high- 
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Figure 2. Top Panel: We show the star formation histories for a subset of 
model galaxies. All histories have been normalized to unity when integrated 
over cosmic time, to make easier the comparison, such a normalization de- 
fines the units of the y-axis. 

Bottom Panel: The BH mass accretion histories for the same model galax- 
ies, normalized as above. The chosen subset of galaxies in our Monte Carlo 
simulation is the same adopted to illustrate the paths in the Mbh ~ ^» plane 
in fig. 1. 



density regions of the primordial density field where the collapse 
and growth of galactic hosts is accelerated. The star formation and 
BH accretion histories of such objects are shown in fig. 2, and can 
be schematically divided in two phases. 

At high redshifts z > 3, the rapid interactions characteriz- 
ing the dense environment of such galaxies are effective in rapidly 
destabilizing the galactic gas and in feeding the central BHs. In fact, 
in such environments and cosmic times, both the interaction rate in 
eq. (2) and the fraction of destabilized gas (eq. 3) are large; the 
first, due to the large densities, and the second due to the large ratio 
m' Im X 1 (i.e., a mass in' of the merging partner close to the mass 
m of the main progenitor) characteristic of this early phase (z > 3) 
when galaxy interactions mainly involve partners with comparable 
mass. The frequent starbursts in star formation and the correspond- 
ing episodes of BH mass accretion at this early cosmic times are 
shown in fig. 2. 

At lower z, the decline of the interaction rate and of the desta- 
bilized fraction /o^^ suppresses the growth of BHs which in our 
model is only due to galaxy interactions (see the drop in the BH 
accretion episodes in the bottom panel of fig. 2), while quiescent 
star formation still proceeds, as is shown by the smooth compo- 
nent in the star formation histories in the top panel of fig. 2. The 
latter, quiescent component of star formation continues to build up 
stellar mass at z < 2 (though at a milder rate), thus lowering the 
Mbh/M, ratio. The overall result is that, when the Mbh/M, ratio is 
normalized to the final local value, the excess T{z) increases with 
redshift. 

Note that the above result is the outcome of two physical pro- 
cesses: 1) the approximatively two-phase growth of cosmic struc- 
tures (with major merging at high redshifts and accretion of small 




Figure 3. The predicted evolution of the BH to stellar mass ratio F (nor- 
malized to the local value, see eq. 1 ), for the evolution of BHs with masses 
Mbh > lO' Mq at z > 4. The 6 filled contours correspond to equally spaced 
values of the fraction of objects with a given value of r{z) at the considered 
redshift: from 0.01 for the lightest filled region to 0.1 for the darkest. For 
reference we draw as a dashed line the local value F(0) = 1 . The data points 
represent the excess F derived from the observed BH masses of the high- 
redshift broad-line AGN following the procedure described in Appendix A. 



clumps at low z) characteristic of hierarchical scenarios, as noted by 
several authors (e.g., Zhao, Jing & Borner 2003; Diemand, Kuhlen 
& Madau 2007; Homan et al. 2007; Ascasibar & Gottloeber 2008); 
2) the interaction-driven scenario for the triggering of AGN activ- 
ity. 

In the above picture, we expect the paths with r(z) > 1 to 
dominate the growth histories of massive objects, formed in biased 
regions of the density field where high-redshift interactions are ex- 
tremely effective, while we expect the high-redshift values of T to 
progressively lower when galaxies formed in less dense environ- 
ment (and hence with a lower mass on average at any given z) are 
considered, due to the progressively lower efficiency of interactions 
in triggering BH accretion at high-z. Indeed, the dependence of the 
above effect on BH mass and on redshift constitutes a typical sig- 
nature of the above interaction-driven model for the growth of BHs 
in galaxies. 



4 RESULTS: THE EVOLUTION OF Mbh - M. FOR 
DIFFERENT AGN POPULATIONS 

4.1 The Early BH Growth and High-redshift QSO 

An efficient way to test the above scenario is to investigate how the 
paths followed by BHs and by their host galaxies to reach the local 
Mbh - M, relation depend on their properties (mass, gas fraction, 
star formation rate), and to compare such specific predictions with 
the corresponding available observations for which we will use the 
mean local value of the Mbh/M, ratio obtained from the Haring & 
Rix (2004) sample to calculate T. We shall devote to this the next 
subsections. 

To explore in detail the growth of BHs in galaxies with differ- 
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Figure 4. Left Panel: The predicted evolution of the BH to stellar mass ratio F (normalized to the local value, see eq. 1), for AGNs selected as to have a 
bolometric luminosity L > 10^-^ erg/s at 1 < z < 2. The 4 filled contours correspond to equally spaced values of the fraction of objects with a given value of 
r{z) at the considered redshift: from 0.01 for the lightest filled region to 0.1 for the darkest. The data points represent the excess T derived from the BH masses 
of the COSMOS broad-line AGN measured by Merloni et al. (2009). Note that the bulk of the BH building up occun'ed ad higher redshift, when most of the 
galaxy interactions triggering the BH accretion take place: the slow evolution T(z) for 1 < z < 2 is driven by the gradual and much milder growth of the stellar 
mass. 

Right Panel: The distribution of log T for all galaxies in the redshift range 1 < z < 2 for model galaxies (solid line) is compared with the corresponding 
observed distribution (Merloni et al. 2009, filled histogram): the latter has been obtained including the lower limits shown in the left panel. The average value 
and the dispersion for log F shown on the top refer to the observed distribution, these values are nearly insensitive to the exclusion of the lower limits. 



ent properties, we start with the most extreme objects, namely, the 
massive BH already in place at high redshifts. 

The time evolution of the Mbh/M, ratio of such extreme ob- 
jects is shown in fig. 3 in terms of the excess T(z) over their local 
value for model BHs with masses Mgjj > 10' Mq at z > 4. For 
such objects r(z) > 1 holds at any time between z = 6 and the 
present, with F » 5 holding at at z = 5-6. These F values are 
consistent with the observations of luminous QSO at 3.6< z <6.4 
(see Appendix A). 

4.2 The BH growth at Intermediate-redshift and BL AGNs 

Now, we focus on BL AGNs at intermediate redshift 1 < z < 2 
and AGN bolometric luminosity L > 10^^ erg/s, for which the 
model results can be compared with detailed existing observations. 
Such objects are expected to form in less biased regions of the pri- 
mordial density field compared to the class of high-redshift QSO 
discussed in the previous section; in fact, they have a comparable 
AGN luminosity but are found at a lower redshifts, and thus are ex- 
pected to originate from merging histories characterized by a lower 
rate of high-redshift encounters. Therefore, according to our pic- 
ture discussed at the end of Sect. 3, we expect such objects to show 
a lower excess F(z). This is indeed the case, as illustrated in fig. 
4, where we compare the model predictions with observations by 
Merloni et al. (2009). 

The distribution of F for such intermediate objects is still dom- 
inated by F > 1 at any refdshift; the model predicts typical values 
Fa; 1.5-2, smaller than those attained by the high-redshift QSO 
shown in fig. 3, and consistent with the observational range. Be- 
sides providing an average F close to the observational value, it is 
interesting to note how the model predicts a scatter in remarkable 



agreement with the observed distribution (see right panel in fig. 
4). This indeed constitute an extremely important test for models 
based on the hierarchical scenarios, since the scatter is directly re- 
lated by the spread in the merging histories of the host galaxies, a 
specific prediction of these models which cannot be tuned though 
adjustable free parameters. 

Note that the model predictions shown in fig. 4 refer to all 
AGNs in the redshift and luminosity ranges specified above, while 
the data we compare with include only unobscured objects. How- 
ever, in our model the obscuration properties of AGNs depend only 
on the amount of gas swept by the blast wave originated by the 
AGN at the time of observation (see Sect. 2.3), i.e., on the de- 
tailed hydrodynamic of the interstellar medium during the last, 
short (At < 5 10' yrs) AGN active phase. Since the F ratios of 
model BHs basically depend on the integrated, past history of the 
simulated BHs and host galaxies, we expect distribution of F for 
obscured and unobscured AGNs to be similar. While we cannot 
perform an exact investigation of this issue (we do not model the 
nuclear obscuration of AGNs), a first-order estimate of the effect of 
selecting only obscured objects can be performed by considering 
in our model only AGNs hosted in galaxies with high gas surface 
densities Eja, > 30 Mq pc"^: even in this case the mean value of 
F in the distribution in fig. 4 is only decreased by a factor « 1.2, 
leaving the scatter unchanged. 

4.3 The Growth of Black Holes in Gas-Rich Galaxies: 
Comparison with Submilhmiter Emitting Galaxies 

In our model, a completely different set of BH accretion histories 
is that corresponding to AGNs in sub-mm galaxies (SMGs). Their 
large gas fractions (/go, > 0.6 relative to the total baryonic mass) 
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Figure 5. The predicted evolution of the BH to stellar mass ratio F (normalized to the local value, see eq. 1), is shown in Panels A and B for AGNs in SMG-like 
galaxies selected from our Monte Carlo simulation according to the criteria i)-iii) (see text). The 5 filled contours correspond to equally spaced values of the 
fraction of objects with a given value of r(^) at the considered redshift: from 0.01 for the lightest filled region to 0.1 for the darkest. Panels A and B refer to 
the different selection in X-ray luminosity of the AGNs shown on the top of the Panels. We compare with data for SMGs with X-ray obscured AGNs (Panel 
A) and with BL AGNs (Panel B) in the same luminosity range adopted for the selection of model galaxies; the data correspond to individual galaxies for 
which both BH and stellar or dynamical masses measurements were available in the literature (see Appendix B). The errorbars in the Panel A show only the 
uncertainties in the data due to the adopted Eddington ratio A, while errorbars on Panel B include the observational uncertainties in both the BH and the stellar 
mass measurements. We also show as a shaded area the 1 - o" uncertainty region around the value of T obtained by Alexander et al. (2008) from average 
BH (adopting A = 0.2) and stellar masses (obtained from near IR luminosity and CO line widths) taken from various samples (see Alexander et al. 2008 and 
references therein). Note that the extremely gas-rich galaxies selected as SMGs represent rare evolutionary paths (see text for the computed and the observed 
number densities) compared to typical galaxies at such redshifts: the statistical fluctuations associated with the low number of selected SMGs in our model 
reflects in the irregular contours and the isolated bins shown in the figure. 



and star formation rate (m, = 100 - 1000 Mg/yr, see Introduction) 
indicate that these galaxies originate from merging histories char- 
acterized by less prominent high-redshft (z > 4) starbursts and BH 
accretion episodes, so that a large fraction of gas is left available at 
lower redshifts z ~ 2 - 3. To compute quantitatively the predicted 
Mbh and M, for such galaxies and to compare with existing ob- 
servations, we selected from our Monte Carlo simulations galaxies 
with i) gas fractions /j^j > 0.7, ii) star formation rates m, > 100 
Mg/yr; iii) AGNs with X-ray luminosities Li < L < L2, where Li 
and L2 have been chosen as to match the selection criteria adopted 
for the observations we compare with (Borys et al. 2005; Alexander 
et al. 2008); the X-ray luminosities (in the band 2-10 keV) have 
been computed from the bolometric luminosities in eq. (4) adopting 
the bolometric correction by Marconi et al. (2004). 

The resulting predicted distribution of the Mg;y/M, ratio are 
shown in figs. 5a and 5b for SMG with AGN luminosity in two 
different ranges (see caption) as to compare with the T estimates 
of six X-ray obscured and three broad line (BL) SMGs for which 
individual estimates of Mgu and M, were available in the litera- 
ture (see the Appendix B for details on the adopted data). Although 
the comparison with observations results are still indicative (due 
to the large uncertainties associated with the Eddington ratios and 
with the measurements of both the BH and the stellar mass, see 
Appendix B), the model predicts for bulk of the SMG population 
values of r < 1, due to the slower growth of BHs at high red- 
shifts. This is due to the particular merging histories corresponding 
to our selection criteria; large residual gas fractions at z ~ 2.5 im- 
ply less frequent high-redshift encounters at z > 4 and less effec- 



tive bursts and BH accretion episodes at these early epochs. Since 
in our model the high-redshift interactions are the only trigger for 
early BH growth, this results in lower values of T(z) at later times 
(z = 2 - 3). Such results do not depend on the details of the se- 
lection criteria i)-ii) that we adopted to extract SMG-like galaxies 
from our Monte Carlo simulations. Indeed, we have verified that 
adopting different thresholds for the gas fractions ( fg„s > 0.6) and 
star formation rates (up to 400 Mg/yr) the fraction of galaxies with 
r <1 is always dominant. 

It must be noted that the observed SMG at z=2 in the right 
panel of fig. 5 has a T value much lower than that predicted by the 
model at the same redshift. In view of the unavoidable uncertainties 
affecting the present observational data (discussed in Appendix B), 
it will be important to confirm the reliability and the statistical sig- 
nificance of such a data point as a marker of the unobscured AGNs 
in SMGs, since it would constitute a severe test for our model which 
predicts a basically similar behaviour of T(z) for obscured and un- 
obscured AGNs. On the other hand, the predicted distribution of F 
for SMG-like galaxies shows a mild dependence with the AGN lu- 
minosity, being larger for luminous AGN. Although the uncertain- 
ties in present data do not allow yet to make precise comparisons, 
this trend constitutes a model prediction which is at least consistent 
with estimates of T based on average values for Mbh and M, in 
SMG galaxies (see the shaded area). 

An interesting feature of our interaction-driven model for BH 
accretion is that the SMGs produced by the model end up in low- 
redshift descendants which lie below the local Mgu — M, relation 
as shown in fig. 6. In other words, the low-redshift ( z < 2) descen- 
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Figure 6. Typical growth paths for BH and stellar masses for our model 
SMG-like galaxies. The dots represent the final location (at z = 0) of the 
descendents of SMG-like galaxies selected according to our criteria i)-iii) 
(see text) in the redshift range 2 < z < 3. The dashed line represent the 
median value of the local Mbh - relation for all model galaxies at z = 0. 

dants of SMGs are predicted to have BHs with low-intermediate 
masses Mbh = - iC Mq (see fig. 6), and the BHs in SMG- 
like galaxies at 2 < z < 3 are in an active growing phase, as it is 
shown in fig. 6 by their steeply rising paths in the M, - Mgu plane. 
Conversely, we computed that a fraction a 10% of local galaxies 
with BH masses in the range Mbh = 10** - 10' Mq had a progenitor 
which passed through an SMG-phase at some redshift in the range 
2<z<3. 

Such a picture is in agreement with the independent findings 
of Alexander et al. (2008), based on the larger number density of 
SMGs compared to that of local galaxies hosting BHs with masses 
exceeding 10' Mq. To check for the full consistency of our results 
and of our interpretation with the observations of SMGs, we have 
compared the number densities of the model SMG galaxies with 
the number densities psmg - 2.5xl0"''Mpc"'' measured by Swin- 
bank et al. (2006) based on the redshift distribution obtained by 
Chapman et al. (2005) for z=l-3.5 SMGs with Sgso,,,,, > 5mJy. To 
perform such a comparison we first converted the above flux limit 
in terms of star formation rate (using the relations in Swinbank et 
al. 2008), since our semi-analytic model does not include the sub- 
mm emission from dust; for the model SMG galaxies above the 
resulting threshold in star formation rate the predicted number den- 
sity in the same redshift interval is Psmg = 4.4xlO"'Mpc"^ when 
SMGs with /ja, > 0.6 are considered, and Psmg = 1.9xlO"^Mpc"^ 
for SMGs with fg^,, > 0.7, fully consistent with the observational 
estimates. 



5 DISCUSSION 

5.1 The growth of BHs for different AGN populations 

The results presented in the previous sections lead to the following 
global picture of our cosmological, interaction-driven model for the 
BH accretion. The relative growth speed of SMBHs compared to 
the stellar mass is characterized by a large spread in the galaxy pop- 
ulation, with an overall trend for a faster BH growth compared to 
galaxies at high redshifts z > 3 ; this is basically due to the fact that 
BH grow only when interactions are effective, i.e., at high redshifts 
z > 3, while star formation proceeds not only through impulsive 
bursts (at high redshifts) but also through quiescent star formation 



which continues to build up stellar mass at z < 2, though at a lower 
rate (see fig. 2a). In terms of the ratio r(z) such a global behaviour 
is represented by fig. 7, where the distribution of T is shown by the 
coloured contours as a function of redshift for the entire (i.e., for all 
galaxies containing BH with masses M > 10^ Mq). A striking fea- 
ture of such a distribution is the modest increase of r(z) for z ^ 3 
as opposed to its rapid upturn at higher redshifts z <; 4; such a mild 
increase of the global T at low-intermediate redshifts is in agree- 
ment with that derived by Merloni et al. (2004) and Hopkins et al. 
(2006) by comparing the observed integrated BH and stellar mass 
densities at such redshifts. 

Note that the overall predicted trend for larger values of T 
with increasing redshift is associated to an increasing the spread 
in the distribution. This constitutes a quantitative prediction for the 
z-dependence of the Lauer et al. (2007) bias. In this respect, the 
r distributions shown in fig. 3, 4, and 5 can be interpreted as a 
quantitative evaluation of the incidence of such a bias in differently 
selected samples at different redshifts. In fact, fig. 7 shows that the 
different observations we compared with in the previous Sections 
are not fully representative of the global behaviour since they re- 
fer to specific populations selected according to different criteria 
which sample only a portion (in the case of SMG actually a minor 
fraction) of the whole galaxy population. 

5.2 Model predictions and downsizing 

A major result of our previous sections is that the evolution of T(z) 
is a strong function of the BH mass. Such a strong mass dependence 
of r(z) is illustrated by fig. 8 (left panel), where we show with the 
contours the average values T(z) as a function of the final (z=0) BH 
mass (represented in the y-axis) and of redshift; the average runs 
over all paths leading to the final BH mass (and hence over all the 
possible main progenitors at a given redshift). The figure illustrates 
that indeed for massive objects Mbh ^ 10' Mq at high redshifts 
z > 4 large values of F > 3 are expected, since such BH form 
in the most biased regions of the density field where high-redshift 
interactions were favoured; values of F > 1.5 - 2 are also natural 
for BH at intermediate redshifts 1 < z < 2, while the low values 
of F observed in SMG at 2 < z < 3 are not representative of the 
whole AGN population at such redhifts (the average F represented 
in fig. 8 is « 1.5), but instead only apply to a population originated 
by peculiar star formation histories (see section 4.3). 

An immediate implication of the above is that massive local 
galaxies and their BHs have formed preferentially through paths 
(in the M, - Mbh plane) passing above the local M, - Mbh re- 
lation; this is illustrated in fig. 8 (right panel) where we show the 
fraction of paths characterized by a given F (in two redshift bins) 
leading to local galaxies with mass M, > 10'^ Mq. The paths with 
F > 1 dominate the statistics not only at high-redshift z = 4 - 5 but 
also -though to a lower extent- at z = 1 - 3. As a consequence, in 
the merging-driven scenario SMGs do not represent typical paths 
leading to local massive galaxies, but rather correspond to peculiar 
paths in the tail of the distribution, namely those selected as to lead 
to gas-rich galaxies at z = 1-3; this interpretation is supported 
by the consistency of the predicted number density of such pecu- 
liar paths with the observed number densities of SMGs, shown and 
discussed at the end of sect. 4.3. 

The above mass dependence of F(z) has straightforward impli- 
cations for the so-called AGN downsizing, i.e., the faster building- 
up of luminous AGNs compared with those with lower luminosities 
which is indicated by observations (see, e.g., Marconi et al. 2004); 
such a downsizing effect is also characterizing the faster drop (2.5 
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Figure 7. The predicted evolution of the BH to stellar mass ratio F for the whole galaxy population in our Monte Carlo simulation; the colour code represents, 
at any redshift z, the fraction of objects with excess Y(z). The average values of F derived from observations of SMGs (Alexander et al. 2008), of BL AGNs 
(Merloni et al. 2009) and high-redshift QSOs (see Maiolino et al. 2007 and Appendix A) are shown as shaded areas covering their observed redshift range and 
observational uncertainties. 
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Figure 8. Left Panel: The average values of F at different redshifts (x-axis) coiresponding to BH with a given final mass (y-axis) are represented as contours. 
The average is taken over all main progenitors in our Monte Carlo merging histories leading to the a given final mass. 

Right Panel: The fraction of paths leading to local massive galaxies with M, > lO'" M© as a function of /ogF is shown for two redshift bins: 4 < z < 5 (shaded 
region), and 1 < z < 3 (solid line). 



dex) from z a 2 to the present of luminous {Lx > 10*' erg/s) AGNs 
compared to the low-luminosity population (Lx < 10*^ erg/s), 
which is observed to follow only a gradual decrease (less than 1 
dex in the number density) from z = 2 to z = (see, e.g., Fiore et 
al. 2003; Hasinger , Miyaji & Schmidt 2005; La Franca et al. 2005). 
In fact, the mass dependence of r(z) represented in fig. 8 implies 
that the massive BH actually grow faster than the low-mass BHs. 



This is illustrated by fig. 9 (left panel), where the contours represent 
the fraction of BH mass formed at a given redshift (in the x-axis), 
for different values of the BH final mass in our model; any given 
fraction of the final mass is assembled at lower redshift for the less 
massive objects with Mgu < 10* Mg with respect to more massive 
objects. We also reproduce as a sequence of points a similar con- 
tour plot produced by Marconi et al. (2004) based on the observed 
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Figure 9. Left Panel: The predicted growth histories of BH with different final mass. For any final BH mass (represented in the y-axis), the 10 filled contours 
correspond to diff'erent values of the average fractional mass Mbh(z)IMbm{z. = 0) of the progenitor BH: Mbh(z)IMbh{i = 0)=0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 
0.7, 0.8, 0.9, from the lightest to the darkest contour We have highlightened with thick black contours the levels corresponding to values 0.9, 0.5, and 0.05, 
to compare with points obtained by Marconi et al. (2004). These were obtained from the analysis of the evolution of the observed AGN luminosity function 
(Ueda et al. 2003), under the assumption of constant Eddington ratio A = \. For such ohsen'ational points, the symbols squares, diamonds and circles, mark 
the contour levels corresponding to the growth of 90 %, 50%, and 5% of the final BH mass, respectively. 

Right Panel: We plot as continuous contours the same growth histories from the analysis by Merloni (2009, see text), assuming that accretion at low m to be 
radiatively inefiicient (RIAF). The dashed contours show the same growth history but assuming for all AGNs the same radiative efficiency (8 %) independently 
of 4. 



mass and luminosity evolution of AGNs, under the assumption of 
fixed Eddington ratio A= \, and that galaxy merging does not affect 
the growth of SMBHs; although both assumptions not hold in our 
model, it is interesting to note that the representation of downsizing 
naturally arising from the interaction-driven model resembles that 
extrapolated from observations. 

The general observed trend of larger assembled fraction for 
increasingly massive BH is quite well accounted for by our 
interaction-driven model, although the model seems to show a 
bending in correspondence of the 50% contour level at z x 1.5 
for masses above 10' Mq which is not present in the data; this 
might constitute a true model inadequacy, or simply result from 
the simple assumptions of constant A adopted in the derivation of 
the data points by Marconi et al. (2004) from the analysis of the 
local BH mass function and of the evolution of the AGN lumi- 
nosity functions. To clarify this issue we reproduce in fig. 9b the 
results from an improved version (Merloni 2009, in prep.) of the 
Merloni & Heinz (2008) analysis, based on updated AGN luminos- 
ity functions. Such an analysis is similar in methodology to that by 
Marconi et al. (2004), but avoids the assumption /i = 1 by deriving 
phenomenological, physically motivated relations between the ac- 
cretion rate, the BH mass and the AGN luminosity, and solving for 
a continuity equation for the SMBH mass function evolution. Be- 
yond the details concerning the exact values of the assembled mass 
fraction derived from the observed luminosity functions (which is 
still largely affected by uncertainties, as shown by the differences 
among the results obtained under different assumptions, see cap- 
tion to fig. 9), the comparison with the Merloni & Heinz (2008) 
approach shows i) that the BH growth inferred from the evolu- 
tion of the AGN luminosity functions generally yields a stronger 
downsizing effect at large masses compared with our model, and 
ii) that such an effect is probably not due to the assumptions on A 
adopted to derive the BH growth histories from the observed lumi- 



nosity functions. Thus, the discrepancy is likely to be due to the 
assumption (common to all the approaches based on the observed 
luminosity functions) that the BH growth is always dominated by 
accretion, and that the BH growth due to merging is negligible at 
any z- 

This conclusion is not unexpected when one considers the ba- 
sic features of the hierarchical clustering picture. Indeed, in all ap- 
proached based on hierarchical scenarios, the merging sector of the 
models predicts massive objects to form later; such a trend is in- 
verted when one consider how and when baryons are converted into 
stars or accreted into SMBHs. The overall downsizing effect in hi- 
erarchical models thus results from the balance between the two 
effects above; the downsizing trend resulting at low-intermediate 
masses in our model means that accretion largely dominates over 
merging up to very large final BH masses Mbh ^ 10' M©. 

While the comparison performed in Fig. 9 is just a first step 
toward a detailed observational testing of the downsizing proper- 
ties of cosmological models for the growth of SMBHs, the global 
picture emerging from the interaction-driven model for the growth 
of SMBHs seems to consistently account for a wide set of indepen- 
dent observations concerning the growth of SMBHs; these range 
from the distant QSOs, to the intermediate-z BL AGNs, to the ex- 
treme population constituted by the SMGs at z ~2, to the average 
growth histories estimated from the observed luminosity functions, 
extending over a wide range of redshifts < z < 4 and of BH 
masses 10' < Mbh/Mq < 10'. 



6 CONCLUSIONS 

We have investigated the growth of SMBHs relative to the stel- 
lar content of the host galaxy in the framework of an interaction- 
driven model for the feeding of BHs during the merging history 
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of the host galaxies. Both the star formation and the BH accre- 
tion (i.e., the conversion of baryons from galactic gas into stars or 
into BHs) are ultimately determined by the histories of the host 
galaxy potential wells, strongly dependent on their total DM mass 
as predicted by galaxy formation models in a cosmological con- 
text. To test such a model against observations, we worked out spe- 
cific predictions for sub samples of the simulated galaxy population 
corresponding to different observational samples for which obser- 
vational estimates of black hole masses Mbh and stellar masses 
M, were available. Specifically, for the evolution of the ratio T = 
(Mbh/M,)(z)/(Mbh/M,)(z = 0) we found that: 

• Massive local galaxies (M, > 10'^ Mq) and their BHs have 
formed preferentially through paths (in the M, — Mbh plane) pass- 
ing above the local M» - Mbh relation. 

• The average value and the spread of T(z) increase over the local 
value for increasing redshift (see fig. 7) . Such an increase is larger 
for massive BH (fig. 8), originating from biased overdense regions 
of the primordial density field. This constitutes a specific prediction 
of the interaction-driven models for the growth of SMBHs, since it 
is due to the effectiveness of interactions in triggering BH accre- 
tion in high-density environments at high redshifts z > 4. In our 
model, the galaxy interactions constitute the only trigger for the BH 
growth, while star formation can proceed in a quiescent model at 
lower redshifts even in the absence of starbursts. Therefore, at high 
redshifts interactions rapidly trigger the accretion of BHs and the 
impulsive star formation, while at lower redshifts z g 2 the interac- 
tion rate drops, and only star formation can proceed in the quiescent 
model, leading to a decrease of T. The larger is the galaxy and BH 
mass, the faster is the gas conversion into stars and the BH growth 
at high redshift, the stronger the decUne of T at low z. 

• Due to the above physical mechanism, extremely massive BHs 
in high redshift QSOs at z > 4 are predicted to have typical val- 
ues of r ~ 4 (see fig. 3), consistently with observational estimates 
(e.g., Willott et al. 2003; McLure & Dunlop 2004; Vestergaard et al. 
2004, Walter et al. 2004; Maiolino et al. 2007; Riechers et al. 2008); 
this corresponds to BHs growing faster than their host galaxy stel- 
lar mass in massive galaxies at high redshifts. At lower z, lower 
mass BH hosted in luminous BL AGNs are predicted to have lower 
values of r « 2 (see fig.4). A direct implication of this is that AGN 
feedback begins to affect the galaxy properies already at high red- 
shifts z > 4. 

• In such a framework, SMG galaxies correspond to the rare merg- 
ing histories leading to i) a large fraction of gas available for impul- 
sive star formation at z « 2, ii) major merging events at such red- 
shifts. These peculiar merging histories are characterized by lower 
interaction rates at high redshifts, which otherwise would have ex- 
hausted the galactic cold gas reservoirs; in turn, this corresponds 
to lower values of F « 0.3 - 1 for such objects (fig. 5), which are 
predicted to be building up a substantial fraction of their final BH 
mass at such redshifts z ~ 2 - 3. In addition, this objects are pre- 
dicted to end up their evolution as low-to-intermediate mass BHs 
with Mbh :£ 10' M^ (fig. 6); such a prediction is in agreement with 
the independent findings of Alexander et al. (2008), based on the 
large number density of SMGs which is larger than that of local 
galaxies hosting BHs with masses exceeding 10' M^; in fact, the 
observed SMG density is consistent with the model predictions, as 
shown in Sect. 4.3. 

The above comparison with observations probe the model in 
different range of masses and redshifts. The global picture emerg- 
ing from the model (see Discussion in Sect. 5) is consistent with a 
downsizing scenario, where massive BHs accrete a larger fraction 
of their final mass at high redshifts z > 4; this is a natural out- 



come of cosmological interaction-driven model for the growth of 
BHs, due to the larger effectiveness of interactions in high-density 
enviroimients, which constitute the density regions where massive 
objects form. Note that this is entirely consistent with the hierar- 
chical building up of cosmic structures, since the number density 
of such massive objects is however increasing with time. 
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APPENDIX A: THE OBSERVED BLACK HOLE-GALAXY 
MASS RATIO AT 3.9< Z <6.4 

A detailed discussion on the BH and galaxy masses in QSO at z 
>4 will be presented in a forthcoming paper by Riechers et al. 
(priv. comm.). In this appendix we obtain estimates of the BH 
and galaxy masses for a sample of four QSO at z <;4 based on 
informations already available in the literature. The T values at 
3.9<z<6.4 in Figs.l, 3 and 7, were obtained from a few quasars 
whose dynamical mass of the host galaxy could be inferred by 
high angular resolution CO milhmctric observations. In paricular 
the CO kinematical data were obtained from Walter et al. (2004), 
for SDSSl 148+5251 at z=6.42, from Maiolino et al. (2007), 
for SDSS0338+0021 at z=5.03, from Riechers et al. (2008), for 
BRn335-0417 at z=4.407, and from Riechers et al. (2009), for 
APM08279+5255 at z=3.911. Note that for SDSS0338+0021 
the dynamical mass is actually only an upper limit, since the CO 
emission is not resolved. We corrected the inferred dynamical 
masses for an average inclination of the molecular disc of 15° in 
agreement with the bias inferred by Carilli & Wang (2006) for 
the inclination of the discs in quasar hosts. We inferred the stellar 
mass by subtracting from the dynamical mass the gas mass (as 
reported in the same references reporting the CO observations) 
and the black hole mass discussed below. Black hole masses 
were inferred for the same quasars by exploiting virial estimators 
based on the width of the broad emission lines and the continuum 
luminosity, calibrated on local AGNs (e.g. Vestergaard & Peterson 
2006, Onken et al. 2004). In particular black hole masses were 
estimated (by using optical and near-lR spectra) by the following 
authors: Barth, Martini & Nelson (2003), for SDSSl 148+5251, 
Dietrich & Hamann (2004) for SDSS0338+0021 (note that the BH 
mass in Maiolino et al. 2007 was underestimated due to calibration 
problems). Shields et al. (2006) for RI1335-0417, Riechers et 
al. (2009) for APM08279+5255. For what concerns the lensing 
factor of APM08279+5255, we adopted the value yu=4 for the CO 
emission in the host galaxy, as inferred by Riechers et al. (2009). 
As discussed by the same authors, the lensing factor of the nuclear 
emission, which is important to estimate the BH mass, can be 
significantly larger, if the latter is close to the caustic. To calculate 
the BH mass of this source we take ju„„c ~ 20, which is the the 
logarithmic average of the former value and the previous estimates 
based on the analysis of the nuclear emission lensed image (Ibata 
et al. 1999, Egami et al. 2000, Munoz et al. 2001, Lewis et al. 2002). 

The BH masses were corrected for the effect of radiation pres- 
sure on the viral estimates, which is particularly important for such 
powerful QSOs, by following the prescription given in Marconi et 
al. (2008). 



Greene & Ho (2005) virial black hole mass estimator. We have as- 
signed an uncertainty of 0.5 dex to these estimates, equal to the Icr 
uncertainty of the virial relation adopted by the authors. 
The BH masses of the X-ray obscured SMG are derived from 
Alexander et al. (2005). They estimated Mbh by converting the X- 
ray luminosity in the (0.5-8) keV band to the AGN bolometric lu- 
minosity (assuming a bolometric correction of 6*^^ %) under the 
assumption that the accretion is Eddington limited. In fig. 5a the 
vertical bars correspond to the Mbh we obtain assuming Eddington 
ratio (A=hi,„ifLEjj) of 0.1, 0.2 and 1 . 

The stellar masses, M,, of the obscured SMGs and of one BL 
SMG (SMGJ123635+621424, star in fig. 5b ) are from Borys et al. 
(2005). They estimated stellar masses using rest-frame UV/near-IR 
luminosities and a mass-to-light ratio of 3.2 for these stellar popu- 
lation. The uncertainties in these M» estimates are from 0.01 dex to 
0.23 dex (see Borys et al. 2005). 

The r values for the other two BL SMG 
(SMGJ131222.3+423814 and SMGJ163706.5+405313, dia- 
monds in fig. 5b) are obtained using dynamical mass estimates, 
which are derived from CO line width measurements (Coppin 
et al. 2008; Greve et al. 2005). CO line widths can be converted 
into dynamical masses assuming a size and inclination for the gas 
reservoir. These estimates assumed a spatial extent of the CO gas 
of R- 2 Kpc. The main uncertainties of these dynamical mass 
estimates are due to the unknown extention of the CO region and 
of the unknown inclination angle. It should also to be noted that 
the dynamical mass includes gas and stellar masses, therefore the 
diamonds in fig. 5b represent lower limits for F. 



APPENDIX B: THE BLACK HOLE AND STELLAR MASS 
MEASUREMENTS IN INDIVIDUAL SMG GALAXIES 

In figs. 5a and 5b the F values of six X-ray obscured and three 
BL SMG galaxies are reported. The BH mass estimates of the BL 
SMGs are taken from Alexander et al. (2008). They estimated Mbh 
from the luminosity and line width of the broad Hor line and the 



